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Abstract. A new approximation hierarchy, called the LPSUBm scheme, is described 
for the coupled cluster method (CCM). It is applicable to systems defined on a regular 
spatial lattice. We then apply it to two well-studied prototypical (spin-1/2 Heisenberg 
antiferromagnetic) spin-lattice models, namely: the XX Z and the XY models on the 
square lattice in two dimensions. Results are obtained in each case for the ground- 
state energy, the ground-state sublattice magnetization and the quantum critical point. 
They are all in good agreement with those from such alternative methods as spin- wave 
theory, series expansions, quantum Monte Carlo methods and the CCM using the 
alternative LSUBm and DSUBm schemes. Each of the three CCM schemes (LSUBm, 
DSUBm and LPSUBm) for use with systems defined on a regular spatial lattice is 
shown to have its own advantages in particular applications. 

PACS numbers: 75.10.Jm, 75.30.Gw, 75.30.Kz, 75.50.Ee 
1. Introduction 

The coupled cluster method (CCM) [TH9] is widely recognized nowadays as providing 
one of the most powerful, most universally applicable, and numerically most accurate 
at attainable levels of computational implementation, of all available ab initio methods 
of microscopic quantum many-body theory. The number of successful applications of 
the CCM to a wide range of physical and chemical systems is now impressively large. 
Some typical examples, from among many others, of systems existing in the spatial 
continuum, and to which the method has been applied, include the electron gas [T0HT2] . 
atomic nuclei and nuclear matter [13l[TI], and molecules |15] . In these and many other 
cases the numerical results obtained with the CCM are either the best or among the 
best available. For the case of the electron gas, for example, which is still one of the 
most intensely studied of all quantum many-body systems, the CCM results [12] for the 
correlation energy agree over the entire metallic density range to within less than one 
millihartree per electron (i.e., better than 1%) with the essentially exact Green's function 
Monte Carlo results available for this system [16]. More recently and more relevantly 
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for the present discussion, the CCM has also been very successfully applied to systems 
on a discrete spatial lattice, such as spin-lattice models of quantum magnetism [TTl - f2"9] . 

One of the features of the CCM, in which it differs from many other techniques 
for dealing with quantum many-body systems, is that, if required, it deals from the 
outset with infinite systems. Thus, one never needs to take explicitly the limit iV — > oo, 
where N is the number of interacting particles or the number of lattice sites. On 
the other hand, of course, the method does require us to make approximations for its 
implementation. These typically involve making selections for which terms to include 
in the cluster expansions for the correlation operators that are intrinsic to the way the 
method parametrizes the many-body wave functions, as we describe more fully in section 
|2] below. 

We and our collaborators have developed previously several efficient and systematic 
approximation schemes for the CCM that are specifically geared to use with lattice 
systems [T8 | I30H33] . The most widely used and the most successful such CCM 
approximation schemes for spin-lattice systems up to now have been the so-called 
LSUBm, and SUBn-m. schemes discussed in detail below in section HI The LSUBm 
scheme in particular has been demonstrated on many occasions to be highly accurate 
in practice for a wide variety of strongly correlated spin systems. Of special importance 
is the fact that the scheme seems to be equally applicable to both frustrated and 
unfrustrated systems, with comparable levels of accuracy attained in both cases. 
Nevertheless, a disadvantage of the LSUBm scheme is that the number of spin 
configurations retained at a given level in describing the many-body correlations present 
in the wave functions, rises very rapidly (and typically super-exponentially) with the 
truncation index m. Since we typically then have to take the limit m — > oo numerically 
to obtain estimates for exact physical properties of the system, it is desirable to have 
calculations at as many values of the truncation index m as possible. 

This one drawback of the prevailing LSUBm scheme has led us recently to develop 
an alternative scheme, the so-called DSUBm scheme [33]. A primary aim of any such 
new scheme should be that in practical applications of it one is able to implement 
more levels of approximation (i.e., to use more values of the index m) than in the 
corresponding LSUBm scheme for the same problem. In this way one thus has more 
data points available for the necessary m — > oo extrapolations, for calculated physical 
quantities, to the exact limit where all spin configurations are retained in the many-body 
wave functions. A second very desirable feature of any such new scheme is that it also 
captures the physically most important multi-spin configurations in the system wave 
functions at relatively low orders in the index m, so that physical properties converge 
more rapidly as m is increased. 

Although the recently developed DSUBm scheme [33] partially met the above 
criteria, there is no doubt that users of the CCM would still welcome more choices 
of approximation schemes. In that context the principal aim of the present work is to 
outline a further such scheme that we now denote as the LPSUBm scheme, and which is 
also specifically designed for use with lattice systems. The scheme is both motivated on 
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physical grounds and its merits illustrated by applications to some stereotypical models 
that have been well studied previously by other techniques, including the CCM itself 
but with other approximation schemes. 

The general formalism of the CCM is first briefly outlined in section |2j after which 
we discuss its specific applications to systems confined to the sites of a regular spatial 
lattice in section [3j In section H] we first describe the existing CCM truncation schemes 
for spin-lattice systems, and then motivate and describe the new LPSUBm scheme. 
The accuracy of the new scheme in practice is then illustrated by applying it to two 
well-studied antiferromagnetic spin-lattice models J20j[22] , namely the spin-half XX Z 
and XY models on the two-dimensional (2D) square lattice. Both models contain a 
free parameter in the Hamiltonian which, as it is varied, carries the zero-temperature 
models through a quantum phase transition at some critical value of this parameter. 
Both models have previously been the subject of CCM studies, using the LSUBm and 
DSUBm truncation schemes, to calculate the ground-state (gs) energy and gs order 
parameter (which, in the present cases, is the sublattice magnetization). 

We note that all microscopic techniques applied to infinite spin-lattice problems 
need to be extrapolated in terms of some appropriate parameter. For example, for such 
main alternative methods to the CCM as the exact diagonalization of small clusters and 
quantum Monte Carlo simulations of larger clusters, the extrapolation parameter is the 
number of lattice sites N. As previously noted, one huge advantage of the CCM is that 
it exactly preserves the Goldstone linked-cluster theorem, and hence size extensively, at 
all levels of approximation. Hence we may (and do) work in the limit of infinite lattice 
size (N — > oo) from the very beginning. By contrast, the extrapolations for the CCM 
are done in terms of some truncation index m, where in the limit m — > oo we retain all 
possible spin configurations in the wave functions of the system, and the calculations 
become formally exact. The extrapolation schemes used in practice [22 | l3"T | l3"3]-f36J are 
themselves also first described in section [5j The new LPSUBm scheme is then applied 
to the spin-half XXZ model and the spin-half XY model, both on the 2D square 
lattice, in sections |6] and [7] respectively. Results are compared critically with those 
from corresponding CCM studies using the alternative LSUBm and DSUBm schemes, 
as well as with the best results from other methods. We conclude in section [8] with a 
summary and discussion of our main findings. 

2. Review of the CCM formalism 

We first briefly describe the CCM formalism. The interested reader is referred, for 
example, to Refs. [8,9J for further details. In any application of the CCM a first step 
is to choose a normalized model (or reference) state |$) that can act as a cyclic vector 
with respect to a complete set of mutually commuting multi-configurational creation 
operators Cf = (Cjy. The index / here is a set-index that labels and uniquely identifies 
the many-particle configuration created in the state Cf |$). The exact ket and bra gs 
energy eigenstates and of the many-body system are then parametrized in the 



The LPSUBm Approximation Scheme for the Coupled Cluster Method 



4 



CCM form as: 

|*>=e 5 |$>; S = £S,C+ (1) 



(^| = ($|5e- 5 ; S=1+^2SjCj, (2) 



where 



H\V) = E\V); {V\H = E(% (3) 

where we have defined Cq = 1 = Cq . The requirements on the multi-configurational 
creation operators are that any many-particle state can be written exactly and uniquely 
as a linear combination of the states {C/|$)}, which hence fulfill the completeness 
relation 

52C+\Q)(Q\Cf = 1 = |$}($| +£<7+|*)<*|C7, (4) 

I IjiO 

together with the conditions, 

C / -|$) = 0=($|C+; (5) 

[Ct,C+] = = [CT,Cj]. (6) 

In practice approximations are necessary to restrict the label set / to some finite 
(e.g., LSUBm) or infinite (e.g., SUBn) subset, as described more fully below. The 
correlation operator S is a linked-cluster operator and is decomposed in terms of a 
complete set of creation operators Cf. When acting on the model state it creates 
excitations that are correlated cluster states. Although the manifest Hermiticity, 
((^|) f = |\&)/(^|^), is lost, the normalization conditions (*|*) = ($|^) = ($1$) = 1 
are preserved. The CCM Schrodinger equations ([3]) are thus written as 

#e 5 |$) = £e 5 |$); {<$>\Se~ s H = E($\Se- s ; (7) 

and their equivalent similarity-transformed forms become 

e- s He s \$) = E\$); (^\Se~ s He s = E($\S. (8) 

While the parametrizations of equations ([I]) and (ED are not manifestly Hermitian 
conjugate, it is very important to note that they do preserve the important Hellmann- 
Feynman theorem at any level of approximations (viz., under any truncation of the 
complete set of many-particle configurations I) [9]. Furthermore, the amplitudes (Si, Si) 
form canonically conjugate pairs in a time- dependent version of the CCM, by contrast 
with the pairs (Si, Sj) coming from a manifestly Hermitian-conjugate representation for 

= (($|e 5t e 5 |$))~ 1 ($|e 5t , which are not canonically conjugate to one another [6]. 

The static gs CCM correlation operators S and S contain the real c-number 
correlation coefficients Si and Si that need to be calculated. Clearly, once the coefficients 
{Si, Si} are known, all other gs properties of the many-body system can be derived from 
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them. Thus, the gs expectation value of an arbitrary operator A, for example, can be 
expressed as 

A = (A) = = ($|5e- 5 Ae 5 |$) ee A(Si,Si). (9) 

The gs correlation coefficients {iS/,iS/} are now found by simply inserting the 
parametrizations of equations ([1]) and ([2]) into the similarity-transformed Schrodinger 
equations (JHJ), and projecting onto the complete sets of states {($|C7} and {(7/1$)}, 
respectively, 

($|C7e _s #e s |$) =0; VI ^ 0. (10) 

($\S(e- s He s - E)C+\$) = 0; V/ ^ 0. (11) 

By pre- multiplying the ket-state equation (jSj) with the state (^IS'C/ and using the 
commutation relation ([6]) it is easy to show that equation ( Till may be rewritten, in the 
form 

(<S>\Se- s [H, C+]e 5 |$) = 0; V/ ^ 0. (12) 

Equations ( |T0|) -(fT2l) may be equivalently derived by requiring that the gs energy 
expectation value, H ee (ty\H\^) = (<&\Se~ s He s \<&), is minimized with respect to the 
entire set {Si, Si}. In practice we thus need to solve equations (flOl and (fl2"|) for the set 
{Si, Si}. We note that equations (EJ) and (JTOl show that the gs energy at the stationary 
point has the simple form 

E ee £(«S/) = ($|e _s i7e s |$), (13) 

which also follows immediately from the ket-state equation (jSJ) by projecting it onto the 
state ($|. It is important to note, however, that this (bi-) variational formulation does 
not necessarily lead to an upper bound for E when the summations over the index set 
{/} for S and S in equations ([I]) and (j2J) are truncated, due to the lack of manifest 
Hermiticity when such approximations are made. Nevertheless, as we have pointed out 
above, one can prove [9] that the important Hellmann-Feynman theorem is preserved 
in all such approximations. 

Equation (TTUT) now represents a coupled set of multinomial equations for the c- 
number correlation coefficients {Si}. The well-known nested commutator expansion of 
the similarity-transformed Hamiltonian, 

e- s He s = H+[H,S} + S], S] + ■ ■ ■ , (14) 

and the fact that all of the individual components of S in the decomposition of equation 
([T]) commute with one another by construction [and see equation ([6}], together imply 
that each element of S in equation (CQ) is linked directly to the Hamiltonian in each of 
the terms in equation ffT4"|) . Thus, each of the coupled equations ffTUj) is of Goldstone 
linked- cluster type, thereby also guaranteeing that all extensive variables, such as the 
energy, scale linearly with particle number N. Thus, at any level of approximation 
obtained by truncation in the summations on the index / in the parametrizations of 
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equations (pQ) and we may (and, in practice, do) work from the outset in the limit 
N — > oo of an infinite system. 

It is now also important for practical applications to note that each of the seemingly 
infinite-order (in S) linked-cluster equations (11 01) will actually be of finite length when 
expanded using equation (114j) . The reason for this is that the otherwise infinite series in 
equation (j!4p will actually terminate at a finite order, provided only (as is usually 
the case, including those for the Hamiltonians considered in this paper) that each 
term in the Hamiltonian H contains a finite number of single-particle destruction 
operators defined with respect to the reference (or generalized vacuum) state |$). In 
this way the CCM parametrization naturally leads to a workable scheme, that can be 
implemented computationally in an efficient manner, to evaluate the set of configuration 
coefficients {Si, Si} by solving the coupled sets of equations (fTUl) and ( fl2j) . once we have 
devised practical and systematic truncation hierarchies for limiting the set of multi- 
configurational set-indices {/} to some suitable finite or infinite subset. We turn our 
attention to such truncation schemes in section H] after first reviewing the application 
of the method (described in general terms above) to the specific case of spin-lattice 
systems. 

3. Review of the CCM for spin-lattice systems 

We now briefly describe how the general CCM formalism outlined in section [2] is 
implemented for spin-lattice problems in practice. As we have already asserted is the 
case for any application of the CCM to a general quantum many-body system, a first 
step is to choose a suitable reference state |$) in which the the state of the spin (viz., in 
practice, its projection onto a specific quantization axis in spin space) on every lattice 
site k is characterized. The choice of 1$) will clearly depend on both the system being 
studied and, more importantly, which of its possible phases is being considered. We 
describe examples of such choices later for the particular models that we utilize here as 
test cases for our new truncation scheme. 

We note firstly that, whatever choice for |<3>) is made, it is very convenient, to 
treat the spins on every lattice site in an arbitrarily given model state |$) as being 
equivalent, in order to create as universal a methodology as possible. A suitably simple 
way of doing so is to introduce a different local quantization axis and a correspondingly 
different set of spin coordinates on each lattice site k, so that all spins, whatever their 
original orientation in |<I>) in the original global spin-coordinate system, align along the 
same direction (which, in order to be definite, we henceforth choose as the negative z 
direction) in these local spin-coordinate frames. In practice this can always be done by 
defining a suitable rotation in spin space of the global spin coordinates at each lattice site 
k. Such rotations are canonical transformations that leave unchanged the fundamental 
spin commutation relations, 



S fc = S k ^ ^ S fc' 



(15) 
(16) 
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among the usual SU(2) spin operators s k = {s%, Sf., sf.) on lattice site k. Each spin has 
a total spin quantum number, where = Sk{sk + 1) is the SU(2) Casimir operator. 
For the models considered here, s/. = s = 1/2, at every lattice site k. 

It is clear that after the local spin axes have been chosen as described above, the 
model state thus has all spins pointing downwards (i.e., in the negative z-direction, 
where z is the quantization axis), 

N 

|$) = | Dk] in the local spin axes, (17) 

k=l 

where | \) = \s, —s) in the usual \s,m s ) notation for single spin states. 

The configuration indices I now simply become a set of lattice site indices, 
I (hi, k 2 , ■ ■ ■ , k m ), and in the local spin frames defined above the corresponding 
generalized multi-configurational creation operators Cf thus become simple products 
of single spin-raising operators, Cf — > s^s^ • • • st ■ Thus, for example, the ket-state 
CCM correlation operator is expressed as 

N 

S = Y1 E s ^k 2 -k m st 1 st 2 ---si m , (is) 

m=l k\k-2.---k m 

and S is similarly defined in terms of the spin-lowering operators s^. Since the operator 
S acts on the state |<3>), in which all spins point along the negative z-axis in the local 
spin-coordinate frames, every lattice site k{ in equation ([TBI can be repeated up to no 
more than 2s times in each term where it is allowed, since a spin s has only (2s + 1) 
possible projections along the quantization axis. 

The allowed configurations are often further constrained in practical applications 
by symmetries in the problem and by conservation laws. An example of the latter is 
provided by the XX Z model considered below in section [6l for which we can easily 
show that the total z-component of spin, = J2k=i s h i n the original global spin 
coordinates, is a good quantum number since [s^, H] — in this case. Finally, for the 
quasiclassical magnetically ordered states that we calculate here for the models in both 
sections [6] and the order parameter is the sublattice magnetization, M, which is given 
within the local spin coordinates defined above as 

N N 

M = ~n^\ E s ^*> = ~n E<*i 5e " s *» e5 i*>- ( 19 ) 

k=l k=l 

The similarity-transformed Hamiltonian H = e~ s He~ s , and all of the 
corresponding matrix elements in equations (I51)- ([TB1) and equation (IT51) . for example, 
may then be evaluated in the local spin coordinate frames by using the nested 
commutator expansion of equation (THj) . the commutator relations of equation fTTS]) . 
and the simple universal relations 

s fc |$) = 0; Vfc, (20) 

s*l^> = ~|^>; v*, (21) 

that hold at all lattice sites in the local spin frames. 
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4. CCM approximation schemes 

When all many-body configurations I are included in the S and S operators in equations 
(0Q) and ([2]) the CCM formalism is exact. In practice, however, it is necessary to use 
approximation schemes to truncate the correlation operators. The main approximation 
scheme used to date for continuous systems is the so-called SUBn scheme described 
below. For systems defined on a regular periodic spatial lattice, we have a further set 
of approximation schemes which are based on the discrete nature of the lattice, such as 
the SUBn-m, LSUBn and DSUBm schemes described below. The various schemes and 
their definitions for spin-lattice systems are: 

(i) the SUBn scheme, in which only the correlations involving n or fewer spin- 
raising operators for 5* are retained, but with no further restrictions on the spatial 
separations of the spins involved in the configurations; 

(ii) the SUBn-m scheme which includes only the subset of all n-spin-flip configurations 
in the SUBn scheme that are defined over all lattice animals of size < m, where 
a lattice animal is defined as a set of contiguous lattice sites, each of which is 
nearest-neighbour to at least one other in the set; and 

(iii) the LSUBm scheme, which includes all possible multi-spin-flip configurations 
defined over all lattice animals of size < m. The LSUBm scheme is thus equivalent 
to the SUBn-m scheme with n = 2sm, for particles of spin quantum number s. For 
example, for spin-1/2 systems, for which no more than one spin-raising operator, 

can be applied at each site k, LSUBm = SUBm-m. 

(iv) the DSUBm scheme, which is defined to include in the correlation operator S all 
possible configurations of spins involving spin-raising operators where the maximum 
length or distance of any two spins apart is defined by L m , where L m is a vector 
joining sites on the lattice and the index m labels lattice vectors in order of size. 
Hence DSUB1 includes only nearest-neighbour pairs, etc. 

We now turn our attention to the new LPSUBm scheme that uses real paths on the 
lattice to determine the fundamental spin configurations. For the LPSUBm scheme, we 
measure distances, P m , along the sides of the lattice, rather than the distance L m used 
in the DSUBm scheme. For example, for a square lattice, we restrict the size of the 
square-lattice plaquette (i.e., the size of the array) by the longest path (P m ) between 
particles in the array, 

P m = k + l; m = k + l, (22) 

where k and I are the sides of the lattice plaquette in the x and y directions. 

Table [1] illustrates the formulation of the spin-array configurations retained in the 
LPSUBm scheme at the m th level of approximation for a 2D square lattice. Similar tables 
can be constructed for an arbitrary regular lattice in any number of dimensions. It shows, 
for example, that the LPSUB5 approximation on a 2D square lattice involves all clusters 
of spins (and their associated spin-raising operators) for which the real path distance 
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Table 1. Illustration of the formulation of the spin-array configurations retained in 
the LPSUBm scheme on a square lattice at the m th level of approximation, in terms 
of lattice increments k and / along the two sides of the square lattice. 



LPSUBm 


Size of square-lattice rectangular 


Maximum no. 




plaquette or size of array 


of spins 




k x / (with P = k + 1) 






[Note: Number of spins on plaquette k x I is f k + 1) x (I + 1)} 




LPSUB1 


1 x 


2 


LPSUB2 


LPSUB1 + 1x1 


4 


LPSUB3 


LPSUB2 + 3x0 + 2x1 


6 


LPSUB4 


LPSUB3 + 4x0 + 3x1 + 2x2 


9 


LPSUB5 


LPSUB4 + 5x0 + 4x1 + 3x2 


12 


LPSUB6 


LPSUB5 + 6x0 + 5x1 + 4x2 + 3x3 


16 


LPSUB7 


LPSUB6 + 7x0 + 6x1 + 5x2 + 4x3 


20 



between any two spins is less than or equal to 5 (lattice spacings). Clearly the LPSUBm 
and the DSUBm schemes both order the multi-spin configurations in terms, roughly, of 
their compactness, whereas the LSUBm scheme orders them, roughly, according to the 
overall size of the lattice animals (or polyominoes), defined as the number of contiguous 
lattice sites involved. 



5. CCM extrapolation schemes 

Each of the above truncated approximations clearly becomes exact when all possible 
multi-spin cluster configurations are retained, i.e., in the limit as n — > oo and/or m — > 
oo. We have considerable experience, for example, with the appropriate extrapolations 
for the LSUBm scheme [22 | [3T | [341436] . that shows that the gs energy behaves in the 
large-m limit as a power series in 1/m 2 , whereas the order parameter M behaves as 
a power series in 1/m (at least for relatively unfrustrated systems). Initial experience 
with the new LPSUBm scheme shows that it behaves similarly with the scaling laws 

E/N = a + ai (^-] + a 2 ( , (23) 



m 2 / V m 2 



for the gs energy (E/N), and 



M = b + h (-) + b 2 (-) , (24) 

for the staggered magnetization (M), respectively, as we show in more detail below for 
the two examples of the spin-1/2 XX Z and XY models on the 2D square lattice. 

In order to fit well to any fitting formula that contains n unknown parameters, one 
should always have at least [n + 1) data points for a robust and stable fit, and in all 
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our CCM calculations in practice we try our best to obey this primary edict, in so far 
as it is possible to do so with the available computing power. In so far as is possible 
we also try to avoid using the least approximate data points (e.g., LSUBm, SUBm- 
m, DSUBm points with m < 2) since these low-m data points are rather far from the 
corresponding large-m limits. In the ensuing discussion we refer to this as our secondary 
edict. Nevertheless, we do include such points if it is necessary to do so to preserve our 
above primary edict. In these latter cases, however, we are always careful to do some 
other careful consistency checks on the robustness and accuracy of our results. 

In the next two sections we now illustrate the use and power of the new LPSUBm 
scheme by applying it to two prototypical spin-half models defined on the 2D square 
lattice, namely the XX Z model in section [6] and the XY model in section [7J 

6. The spin-1/2 antiferromagnetic XX Z model on the square lattice 

As an illustration of the use of the LPSUBm scheme we first consider its application to 
the spin-1/2 XX Z model on the infinite square lattice. The Hamiltonian of the XX Z 
model, in global spin coordinates, is written as 

H xxz = 5>*sJ + s\s) + As***}, (25) 

where the sum on runs over all nearest-neighbour pairs of sites on the lattice and 
counts each pair only once. Since the square lattice is bipartite, we consider N to be 
even, so that each sublattice contains ~iV spins, and we consider only the case where 
iV — > oo. The Neel state is the ground state (GS) in the trivial Ising limit A — >• oo, 
and a phase transition occurs at A = 1. Indeed, the classical GS demonstrates perfect 
Neel order in the z-direction for A > 1, and a similar perfectly ordered x-y planar Neel 
phase for — 1 < A < 1. For A < — 1 the classical GS is a ferromagnet. 

The case A = 1 is equivalent to the isotropic Heisenberg model, whereas A = 
is equivalent to the isotropic version of the XY model considered in section [7J below. 
The z component of total spin, sf,, is a good quantum number as it commutes with the 
Hamiltonian of equation ( 1251) . Thus one may readily check that [s^, Hxxz] = 0. Our 
interest here is in those values of A for which the GS is an antiferromagnet. 

The CCM treatment of any spin system is initiated by choosing an appropriate 
model state |<&) (for a particular regime), so that a linear combinations of products of 
spin-raising operators can be applied to this state and all possible spin configurations 
are determined. There is never a unique choice of model state |$). Clearly our choice 
should be guided by any physical insight that we can bring to bear on the system or, more 
specifically, to that particular phase of it that is under consideration. In the absence of 
any other insight into the quantum many-body system it is common to be guided by 
the behaviour of the corresponding classical system (i.e., equivalently, the system when 
the spin quantum number s — » oo). The XX Z model under consideration provides just 
such an illustrative example. Thus, for A > 1 the classical Hamiltonian of equation 
( 1251 on the 2D square lattice (and, indeed, on any bipartite lattice) is minimized by 
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a perfectly antiferromagnetically Neel-ordered state in the spin ^-direction. However, 
the classical gs energy is minimized by a Neel-ordered state with spins pointing along 
any direction in the spin x-y plane (say, along the spin x-direction) for — 1 < A < 1. 
Either of these states could be used as a CCM model state |<3>) and both are likely to 
be of value in different regimes of A appropriate to the particular quantum phases that 
mimic the corresponding classical phases. For present illustrative purposes we restrict 
ourselves to the z-aligned Neel state as our choice for |$), written schematically as 
1^) = I ' ' ' tttt ' ' ')> m the global spin axes, where | f) = If, +§) an d | l) = ||, — |) in 
the usual \s,m s ) notation. Such a state is, clearly, likely to be a good starting-point for 
all A > 1, down to the expected phase transition at A = 1 from a z- aligned Neel phase 
to an x-y planar Neel phase. 

As indicated in section [3] it is now convenient to perform a rotation of the axes for 
the up-pointing spins (i.e., those on the sublattice with spins in the positive z-direction) 
by 180° about the spin y-axis, so that |<3>) takes the form given by equation (fTTj) . Under 
this rotation, the spin operators on the original up sub-lattice are transformed as 

s x ->• -s x , s y -> s y , s z ->• -s z . (26) 

The Hamiltonian of equation (125]) may thus be rewritten in these local spin coordinate 

H xxz = -\ "£144 + *T*J + 2As ts-]- (27) 

(hi) 

As in any application of the CCM to spin-lattice systems, we now include in 
our approximations at any given order only those fundamental configurations that 
are distinct under the point and space group symmetries of both the lattice and the 
Hamiltonian. The number, Nf, of such fundamental configurations at any level of 
approximation may be further restricted whenever additional conservation laws come 
into play. For example, in our present case, the XXZ Hamiltonian of equation (1231) 
commutes with the total uniform magnetization, s^ = Ylk=i m ^he global spin 
coordinates, where k runs over all lattice sites. The GS is known to lie in the SJ, — 
subspace, and hence we exclude configurations with an odd number of spins or with 
unequal numbers of spins on the two equivalent sublattices of the bipartite square lattice. 
We show in figure [H the fundamental configurations that are accordingly allowed for the 
LPSUBm approximations for this spin- 1/2 XXZ model on the 2D square lattice, with 
1 < m < 3. We see, for example, that Nf = 9 at the LPSUB3 level of approximation. 

The LPSUBm approximations can readily be implemented for the present spin- 1/2 
XXZ model on the 2D square lattice for all values m < 6 with reasonably modest 
computing power. By comparison, the LSUBm and DSUBm schemes can both be 
implemented with comparable computing resources for all values m < 9. Numerical 
results for the gs energy per spin and the sublattice magnetization are shown in table [2] 
at the isotropic point A = 1 at various levels of approximation, and corresponding 
results for the same quantities are displayed graphically in figures [2] and [3] as functions 
of the anisotropy parameter A. 
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Figure 1. The fundamental configurations for the LPSUBto scheme with m — {1, 2, 3} 
for the spin-1/2 XX Z model on a square lattice in two dimensions. The filled circles 
mark the relative positions of the sites of the square lattice on which the spins are 
flipped with respect to the model state. The unfilled circles represent unflipped sites. 



We also show in table|2]for the isotropic Heisenberg Hamiltonian (A = 1) the results 
for the gs energy and sublattice magnetization using the (quadratic) extrapolation 
schemes of equations ( 1231) and (J21J) respectively of the LPSUBm data, employing various 
subsets of results. Comparison is also made with corresponding LSUBm [371 [38] an d 
DSUBm [33] extrapolation schemes for the same model. The results are generally 
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O termination point at A c 
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Figure 2. (Color online) CCM results for the ground-state energy per spin, E/N, as a 
function of the anisotropy parameter A, of the spin-1/2 XX Z model on the 2D square 
lattice, using various LPSUBm approximations based on the ^-aligned Neel model 
state. The LPSUBm results with m = {2,4,6} are extrapolated using the quadratic 
fit of equation (1231) and shown as the curve LPSUBoo. Aj = magnetization point of 
inflexion, described in the text. 
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Table 2. The ground-state energy per spin (E/N) and sublattice magnetization (M) 
for the spin-1/2 XX Z model on the 2D square lattice, obtained using the CCM 
LPSUBm approximation scheme with 1 < m < 6 at A = 1. Nf is the number 
of fundamental configurations at a given LPSUBm, LSUBm or DSUBm level of 
approximation. Aj = LPSUBm sublattice magnetization point of inflexion. The 
LPSUBm results for odd values of m, even values of m and the whole series of m 
are extrapolated separately. These results are compared to calculations using third- 
order spin- wave theory (SWT) [39], series expansions (SE) [40], exact diagonalization 
(ED) [41], quantum Monte Carlo (QMC) [42], LSUBoo extrapolations of the CCM 
LSUBm approximations [H1I3Z1I3H] and the DSUBoo extrapolations of the CCM 
DSUBm approximations |33j . 



ivietiloa 


AT 


E/N 


M 




A 


Max. no. 






A = 


1 






of spins 


LPSUB1 


1 


-0.64833 


0.421 






2 


LPSUB2 


2 


-0.65311 


0.410 


0.258 




4 


LPSUB3 


9 


-0.66442 


0.379 




0.579 


6 


LPSUB4 


35 


-0.66565 


0.372 


0.586 




8 


LPSUB5 


265 


-0.66761 


0.358 




0.766 


12 


LPSUB6 


2852 


-0.66807 


0.354 


0.735 




16 


LSUB8 


1287 


-0.66817 


0.352 




0.844 


8 


LSUB10 


29605 


-0.66870 


0.345 






10 






Extrapolation 






Based on 


LPSUBoo 




-0.66953 


0.320 






m = {1, 3, 5} 


LPSUBoo 




-0.67004 


0.308 


1.093 




m = {2,4,6} 


LPSUBoo 




-0.66867 


0.328 






2 < m < 6 


LPSUBoo 




-0.67107 


0.288 






3 < m < 6 


DSUBoo 




-0.67082 


0.308 


1.009 




m = {6,8,10} 


LSUBoo 




-0.67029 


0.304 






n = {3,5,7,9} 


LSUBoo 




-0.66966 


0.310 






n = {4,6,8,10} 


LSUBoo 




-0.66962 


0.308 






n = {6,8,10} 


SWT 




-0.66999 


0.3069 








SE 




-0.66930 


0.307 








ED 




-0.67000 


0.317 








QMC 




-0.669437(5) 


0.3070(3) 









observed to agree very well with each other. Excellent agreement of all the CCM 
extrapolations is also obtained with the results from the best of the alternative methods 
for this model, including third-order spin- wave theory (SWT) [59] . linked-cluster series 
expansion (SE) techniques [40J, the extrapolations to infinite lattice size (N — > oo) from 
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Figure 3. (Color online) CCM results for the ground-state sublattice magnetization, 
M, as a function of the anisotropy parameter A, of the spin-1/2 XX Z model on 
the 2D square lattice, using various LPSUBm approximations based on the z-aligned 
Neel model state. The LPSUBm results with m — {2, 4, 6} are extrapolated using 
the quadratic fit of equation ([24| and shown as the curve LPSUBoo. Aj = point of 
inflexion in the curve, shown by arrows in the figure. 

the exact diagonalization (ED) of small lattices [H], and quantum Monte Carlo (QMC) 
calculations for larger lattices [42J. 

As discussed in section [5] we always prefer to have at least 4 LPSUBm calculations 
with different values of the truncation index m, to fit to the three unknown parameters 
of the quadratic fitting expressions for E/N and M in equations (123]) and (124]) . This 
primary edict is not violated if we extrapolate the LPSUBm data using both odd and 
even values of m in the range 2<m<6or3<m<6. We note, however, that if we 
extrapolate using only the three even values m = {2, 4, 6} or using the three odd values 
m = {1,3,5} then we violate both the primary and secondary edicts discussed above. 
Nevertheless, the extrapolated results using the even set m = {2, 4, 6} are seen to be in 
good agreement with those from the alternative methods shown in table [2] 

It has been observed and well documented in the past (and see, e.g., Ref. [37]) that 
the CCM LSUBm results for this model (and many others) for both the gs energy E and 
the sublattice magnetization M show a distinct period-2 "staggering" effect with index 
m, according to whether m is even or odd. As a consequence the LSUBm data for both E 
and M converge differently for the even-m and the odd-m sequences. This is very similar 
to what is also observed very frequently in perturbation theory in corresponding even 
and odd orders (43]. As a rule, therefore, the LSUBm data are generally extrapolated 
separately for even m and for odd values of m, since the staggering makes extrapolations 
using both odd and even values together rather difficult. We show in figure H] our 
LPSUBm results for the gs energy per spin and the sublattice magnetization plotted 
against 1/m 2 and 1/m, respectively, for the case A = 1. The higher odd and even 
m values taken together clearly cluster well in both cases on straight lines, thereby 
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Figure 4. (Color online) Illustration of the odd-even staggered nature with respect to 
the truncation parameter m of the LPSUBm scheme results for the ground-state energy 
per spin, E/N, and sublattice magnetization, M, for the spin-1/2 antiferromagnetic 
XX Z model on the 2D square lattice, for the isotropic limiting case A = 1. The 
LPSUB?ti data are plotted against 1/m 2 for E/N and against 1/m for M. The results 
clearly justify the heuristic extrapolation schemes of equations (f2"5]l and (fM]) . 

justifying a posteriori our heuristic extrapolation fits of equations ( 123]) and ( 1241) . Just 
as in the LSUBm case a small but definite "odd-even staggering" effect is observed in 
the LPSUBm data for both the energy and the sublattice magnetization, although it is 
less pronounced than for the corresponding DSUBm [33J and LSUBm data [37J for this 
model. 

Before discussing our LPSUBm results further for this model we note that the 
comparable LSUBm solutions actually terminate at a critical value A c = A c (m), which 
itself depends on the truncation index m [32]. Such LSUBm termination points are very 
common for many spin-lattice systems. They have been very well documented and their 
origin is clearly understood (and see, e.g., Ref. [32]). Thus, in all such termination 





. 5 






E vs m" 2 


m={l,3,5} X 


^ Extrapolated CCM 


m={l,3,5} - 


E vs m" 


m={ 2,4,6} + 


Extrapolated CCM 


m={2,4,6} 









/? Mvsm" 1 


m={l,3,5} X 


'/ Extrapolated CCM 


m={l,3,5} - 


Mvsm" 1 


m={ 2,4,6} + 


Extrapolated CCM 


m={2,4,6} 



0.2 0.4 0.6 0.8 1 



The LPSUBm Approximation Scheme for the Coupled Cluster Method 



16 



point always arises due to the solution of the CCM equations becoming complex at this 
point, beyond which there exist two branches of entirely unphysical complex conjugate 
solutions [32]. In the region where the solution reflecting the true physical solution is real 
there actually also exists another (unstable) real solution. However, only the (shown) 
upper branch of these two solutions reflects the true (stable) physical GS, whereas the 
lower branch does not. The physical branch is usually easily identified in practice as the 
one which becomes exact in some known (e.g., perturbative) limit. This physical branch 
then meets the corresponding unphysical branch at the corresponding termination point, 
beyond which no real solutions exist. The LSUBm termination points are themselves 
also reflections of the quantum phase transitions in the real system, and may hence be 
used to estimate the position of the phase boundary [32J. 

We note that when the LPSUBm approximations are applied to the XXZ model, 
only the odd LPSUBm levels with m > 3 terminate in the same way as do the 
corresponding LSUBn approximations, as shown in figures [2] and [3j The LPSUBm 
solutions with even values m = {2,4,6} do not terminate. We have no convincing 
explanation for this difference in behaviour for two apparently similar schemes applied 
to the same model. Nevertheless, it is still possible to use our LPSUBm data to extract 
an estimate for the physical phase transition point at which the ^-aligned Neel phase 
terminates. As has been justified and utilized elsewhere [20], a point of inflexion at 
A = Aj in the sublattice magnetization M as a function of A also indicates the onset of 
an instability in the system. Such inflexion points Aj = Aj(m) occur for the even values 
of the LPSUBm approximations, as indicated in table [2] and figure [31 The LPSUB?n 
approximations are thus expected to be unphysical for A < Aj(m), and we hence 
show the corresponding results for the gs energy per spin in figure [2] only for values 
Aj > Aj(m). Heuristically, we find that the magnetization inflexion points Aj(m) 
scale linearly with (1/m) in the large m limit, and the extrapolated results shown 
in table [2] have been performed with, Aj = c a + ci(l/m) + C2(l/m) 2 , commensurate 
with the corresponding fits in 1/m 2 and 1/m for the gs energy per spin and sublattice 
magnetization of equations (1231) and (1241) . respectively. The extrapolated values from 
both the LPSUBm and DSUBm schemes are in excellent agreement with the expected 
phase transition point at A c = 1 between two quasiclassical Neel-ordered phases aligned 
along the spin z- axis (for A > 1) and in some arbitrary direction in the spin x-y-plane 
(for |A| < 1). 

Although we do not do so here, the x-y planar Neel phase could itself also easily be 
investigated by another CCM LPSUBm series of calculations based on a model state |$) 
with perfect Neel ordering in, say, the x-direction. Nevertheless, from our results so far 
we observe that the LPSUBm scheme has at, least partially, fulfilled the expectations 
placed upon it for the present model. Accordingly, we now apply it to the second test 
model of the spin- 1/2 XY model on the 2D square lattice. 
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7. The spin-1/2 XY model on the square lattice 



The Hamiltonian of the XY model [20] in global spin coordinates, is written as 

H XY = £[(1 + A)s*s* + (1 - AKaJ]; -1 < A < 1, (28) 

where the sum on again runs over all nearest-neighbour pairs of lattice sites and 
counts each pair only once. We again consider the case of spin-1/2 particles on each site 
of an infinite 2D square lattice. 

For the classical model described by equation ( 1281) . it is clear that the GS is a 
Neel state in the x-direction for < A < 1 and a Neel state in the ^-direction for 
— 1 < A < 0. Hence, since we only consider the case < A < 1, we choose as our CCM 
model state |<3>) for the quantum XY model a Neel state aligned along the x-direction, 

written schematically as, |<&) = | ■ • ■ < H > ■ ■ ■), in the global spin axes. Clearly the 

case — 1 < A < is readily obtained from the case < A < 1 by interchange of the x- 
and y-axes. 

As before in section [6] we now perform our usual rotation of the spin axes on each 
lattice site so that |$) takes the form given by equation (jTTJ) in the rotated local spin 
coordinate frames. Thus, for the spins on the sublattice where they point in the negative 
^-direction in the global spin axes (i.e., the left-pointing spins) we perform a rotation 
of the spin axes by +90° about the spin y-axis. Similarly, for the spins on the other 
sublattice where they point in the positive x-direction in the global spin axes (i.e., the 
right-pointing spins) we perform a rotation of the spin axes by —90° about the spin 
2/-axis. Under these rotations the spin operators are transformed as 

s x -)■ s z , s v ->■ s v , s z ->■ -s x , left-pointing spins; (29 a) 
s x -> -s z , s y -> s v , s z -> s x , right-pointing spins. (296) 



H X y = Yl 



(30) 



The Hamiltonian of equation ( 12811 may thus be rewritten in the local spin coordinate 
axes defined above as 

1 + A)s z s z - ~(1 - A)( S + S + + s-sj) + - A {1 - A)(stsj + S r<f 

Exactly as in the previous application, we now have to evaluate the fundamental 
configurations that are retained in the CCM correlation operators S and S at each 
LPSUBm level of approximation. Although the point and space group symmetries 
of the square lattice (common to both the XXZ and XY models considered here) 
and the two Hamiltonians of equations ( 1271) and ( 1301) are identical, the numbers Nf 
of fundamental configurations for a given LPSUBm level are now larger (except for 
the case m — 1) for the XY model than for the XXZ model, since the uniform 
magnetization is no longer a good quantum number for the XY model, [Hxy, s^} ^ 0. 
Nevertheless, we note from the form of equation ( |30|l . in which the spin-raising and 
spin-lowering operators appear only in combinations that either raise or lower the 
number of spin flips by two (viz., the sfslf and s^sj combinations, respectively) or 
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Figure 5. The fundamental configurations for the LPSUBm scheme with m = {1, 2, 3} 
for the spin-1/2 XY model on a square lattice in two dimensions. The filled circles 
mark the relative positions of the sites of the square lattice on which the spins are 
flipped with respect to the model state. The unfilled circles represent unflipped sites. 



leave them unchanged (viz., the sfsj and s~sj combinations), it is only necessary for 
the s?p = GS to consider fundamental configurations that contain an even number of 
spins. Thus, the main difference for the XY model over the XXZ model is that we 
must now also consider fundamental configurations in which we drop the restriction for 
the former case of having an equal number of spins on the two equivalent sublattices 
of the bipartite square lattice that was appropriate for the latter case. We show in 
figure [5] the fundamental configurations that are allowed for the spin-1/2 XY model on 
the square lattice for the LPSUBm approximation with 1 < m < 3, and we invite the 
reader to compare with the corresponding fundamental configurations for the spin-1/2 
XXZ model on the same square lattice shown in figure 1. The corresponding numbers 
Nf of fundamental configurations for the XY model are also shown in table |3] for the 
LPSUBm approximations with m < 6. 

We present results for the spin-1/2 XY model on the square lattice in the CCM 
LPSUBm approximations for all values m < 6, all of which can be easily computed with 
very modest computing power. Comparable computing power enables the corresponding 
LSUBm scheme to be implemented for all m < 8. Numerical results for the gs energy per 
spin and sublattice magnetization are shown in table [3] at the isotropic point at A = 
at various levels of approximation, and corresponding results for the same gs quantities 
are shown graphically in figures [6] and [7] as functions of the anisotropy parameter A. 
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Table 3. The ground-state energy per spin (E/N) and sub-lattice magnetization (M) 
for the spin-1/2 XY model on the 2D square lattice, obtained using the CCM LPSUBm 
approximation scheme with l<m<6atA = 0. Nf is the number of fundamental 
configurations at a given level of LPSUBm, LSUBm or DSUBm approximation. A c = 
LPSUBm termination point. The LPSUBm results for odd values of m, even values of 
m and the whole series of m are extrapolated separately. These results are compared 
to calculations using series expansions (SE) [44] ■ the quantum Monte Carlo (QMC) 
method [45] , LSUBoo extrapolations of the CCM LSUBm approximations [201 and the 
DSUBoo extrapolations of the CCM DSUBm approximations [33 . 



Method 
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E/N 
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Max. no. 
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of spins 


LPSUB1 
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-0.54031 


0.475 
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LPSUB2 
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-0.54548 


0.464 


-0.401 
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LPSUB3 


13 


-0.54747 


0.457 


-0.178 


6 


LPSUB4 


72 


-0.54812 


0.453 


-0.107 
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LPSUB5 


557 


-0.54842 


0.450 


-0.072 


12 


LPSUB6 


7410 


-0.54857 


0.448 
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16 


LSUB6 


131 


-0.54833 


0.451 


-0.073 
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LSUB8 


2793 


-0.54862 


0.447 


-0.04 
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Extrapolation 






Based on 


LPSUBoo 




-0.54894 


0.437 


-0.017 


2 < m < 5 


LPSUBoo 




-0.54897 


0.435 


-0.006 


3 < m < 5 


LPSUBoo 




-0.54893 


0.436 
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2 < m < 6 


LPSUBoo 
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LPSUBoo 




-0.54888 
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4 < m < 6 


LPSUBoo 




-0.54899 
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m = {1, 3, 5} 


LPSUBoo 




-0.54893 


0.436 
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m = {2,4,6} 


LSUBoo 




-0.54892 


0.435 


0.00 


n = {4,6,8} 


DSUBoo 




-0.54950 


0.436 




m = {3,5,7,9} 


DSUBoo 




-0.54923 


0.437 


0.011 


m = {5,7,9} 


SE 




-0.5488 


0.436 


0.0 




QMC 




-0.54882(2) 


0.437(2) 







NOTES : 

a The LPSUB1 approximation does not terminate. 

b The spin-flip configurations for the LPSUB6 approximation are sufficiently complicated 
and large in number that calculations have only been calculated at present for A = 
for the LPSUB6 case. 
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Figure 6. (Color online) CCM results for the ground-state energy per spin, E/N, 
as a function of the anisotropy parameter A, of the spin- 1/2 XY model on the 2D 
square lattice obtained using the LPSUBm approximation based on the Neel state 
aligned along any axis in the x-y plane. The LPSUBm results with m = {1, 3, 5} are 
extrapolated using equation (|2"3"|) to give the curve labelled DSUBoo. 




Figure 7. (Color online) CCM results for the ground-state sublattice magnetization, 
M, as a function of the anisotropy parameter A, of the spin-1/2 XY model on the 
2D square lattice obtained using various LPSUBm approximations based on the Neel 
state aligned along any axis in the x-y plane. The LPSUBm results with m = {1, 3, 5} 
are extrapolated using equation to give the curve labelled DSUBoo. 



We also show in table |3]for the isotropic XY Hamiltonian (A = 0) the results for the 
gs energy and sublattice magnetization using the (quadratic) extrapolation schemes of 
equations (1231 and (j2l|) respectively of the LPSUBm data, employing various subsets of 
our results, just as for the XX Z model considered previously. We also compare in table[3] 
the present results with the corresponding CCM LSUBm [20] and DSUBm [33] results 
for the same model. All of the CCM results are clearly in excellent agreement both with 
one another and with the results of the best of the alternative methods available for this 
model, including the linked-cluster series expansion (SE) technique [44] and a quantum 
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Figure 8. (Color online) Illustration of the odd-even staggered nature with respect 
to the truncation parameter m of the LPSUBm scheme results for the ground-state 
energy per spin, E/N, and sublattice magnetization, M, for the spin-1/2 XY model 
on the 2D square lattice, for the isotropic limiting case A = 0. The LPSUBm data are 
plotted against 1/m 2 for E/N and against 1/m for M. The results clearly justify the 
heuristic extrapolation schemes of equations (|23|) and (|24|) . 



Monte Carlo (QMC) method [45]. 

We show in figure |S]our LPSUBm results for the present XY model for the gs energy 
per spin and the sublattice magnetization, plotted respectively against 1/m 2 and 1/m, 
for the case A = 0. As previously for the XXZ model, the higher m values cluster 
well on straight lines in both cases, thereby justifying once more our heuristic choice of 
extrapolation fits indicated in equations ( 1231) and ( |24|) . Figures 8(a) and 8(b) once more 



show an "odd-even" staggering effect in the termination index m for the LPSUBm data 
and we have again shown separate extrapolations of our LPSUBm results in table |3] for 
the even-m data and the odd-m data, as well as results using all (higher) values of m. 
It is interesting to note, however, that the staggering effect for this XY model is far 
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less pronounced than for the similar XX Z model in section O We have no compelling 
argument to explain this difference. 

It is interesting to note that for the present XY model the CCM LPSUBm solutions 
(with our choice of model state as a Neel state in the x-direction) now do physically 
terminate for all values of the truncation index m > 1 at a critical value A c = A c (m), 
exactly as commonly occurs (as for the present model) for the LSUBm calculations, 
as we explained above in section [6j Why such LPSUBm terminations occur for all 
values m > 1 for the XY model but not for odd values of m for the previous XXZ 
model is not obvious to us. The corresponding termination points, A c = A c (m), at 
various LPSUBm, LSUBm and DSUBm levels of approximation are shown in table El 
It has been shown previously [31] that A c (m) scales well with (1/m) 2 for the LSUBm 
data, and the LSUBoo result [20] shown in table |3] was obtained by the scaling law, 
A c (m) = do + di(l/m) 2 + ^(l/^) 4 - We find heuristically that the best large-m 
asymptotic behaviour of the LPSUBm data for A c (m) is also against (1/m) 2 as the 
scaling parameter. Accordingly, the LPSUBoo values for A c in table |3] are obtained 
with the same (quadratic) fit, A c (m) = d + <ii(l/m) 2 + (i2(l/m) 4 . We see that both the 
LSUBoo and LPSUBoo results for A c = A c (oo) agree very well with the value A c = 
that is known to be the correct value for the phase transition in the one- dimensional spin- 
1/2 XY chain from the known exact solution [16], and which is believed on symmetry 
grounds also to be the phase transition point for higher dimensions, including the present 
2D square lattice. 

8. Conclusions 

From the two nontrivial benchmark spin-lattice problems that we have investigated 
here, it is clear that the new LPSUBm approximation scheme works well for calculating 
their gs properties and phase boundaries. We have utilized here only the simplest 
extrapolation schemes in the pertinent scaling variables, and have shown that these may 
be chosen, for example, as 1/m 2 for the gs energy and 1/m for the order parameter. 
For further use of the scheme for more complex lattice models (e.g., those exhibiting 
geometric or dynamic frustration) it will be necessary to re-visit the validity of these 
expansions, but a great deal of previous experience in such cases for the LSUBm scheme 
should provide good guidance. 

On the basis of the test results presented here, the LPSUBm scheme clearly does 
not fulfill the first of our two main criteria for introducing it, since the number of 
fundamental configurations, Nf, actually increases even more rapidly with truncation 
index m than for the corresponding LSUBm series of approximations. Nevertheless, our 
second criterion of capturing the physically most important configurations at relatively 
low levels of approximation does seem to be fulfilled, according to our experience with 
the convergence of the LPSUBm sequences for observable quantities. At the very least 
we now have three schemes (LSUBm, DSUBm and LPSUBm) available to us for future 
investigations of more complicated spin-lattice models, each of which has its own merits, 
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and which thus allows us more freedom in future applications of the CCM to quantum 
magnetism. 

It is particularly worth noting too that our preliminary calculations have shown 
that the different schemes show markedly varying patterns of odd-even staggering, both 
for a given scheme applied to different models and for different schemes applied to the 
same model. It is difficult to predict in advance how strong or weak the effect will be for 
a given scheme applied to a specific model. Nevertheless, when the effect is weak one 
can confidently extrapolate the results using both odd and even values of the truncation 
parameter m simultaneously, thereby effectively doubling the number of data points for 
the fit. In such cases our first criterion for an improved scheme has effectively been 
realized over one where the staggering effect is much more pronounced, even though 
the number of fundamental configurations, Nf, may indeed increase more rapidly with 
truncation index m for the former ( "improved" ) scheme than for the latter. 
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